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Evaporation induced flow inside circular wells 
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Abstract. Flow field and height averaged radial velocity inside a droplet evaporating in an open circular well were 
calculated for different modes of liquid evaporation. 

PACS. 47.55.nb Capillary and fhermocapillary flows 



1 Introduction 

Wide used technique of patterning surfaces with solid particles 
utilizes the evaporation of colloidal droplets from a substrate. 
In particular, evaporation of liquid samples is a key problem in 
the development of microarray technology (including labs-on-a 
chip), especially in the case of open reactors [1]. 

A plenty of works are devoted to the evaporation of the 
sessile droplets [2-5]. One assumes that a droplet has a shape 
like a spherical cap. Nevertheless, investigation of the drops 
with more complex geometry is of interest. Thus when a high 
concentrated colloidal droplet evaporates from a substrate, the 
solute particles in solution will form a ring like deposit wall 
near the edge of the drop (Fig. Q~|i. One can consider this case 
as evaporation of a liquid drop inside open circular well with 
vertical walls formed by gel. 

In the paper [1] the flow field inside a liquid sample in 
very thin circular wells was measured. The confidence between 
measured velocity [1] and the estimations obtained with theory 
of ring formation [2] is reasonable. 

In this paper, the relationship between hydrodynamical flow 
inside a liquid sample evaporating in open circular well and the 
mode of evaporation is examined. 

We found analytical expressions of height averaged radial 
velocity for different modes of evaporation. Velocity field in- 
side the circular well was found for the particular case of flat 
air-liquid interface. 
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Fig. 1. Time evolution of the deposit phase growth [6] 



2 Velocity field inside circular well 

2.1 Height averaged radial velocity 

We will consider one particular but interesting case, because of 
its practical applications, when a liquid inside a circular well 
may be described as a thin film. So, in paper [1] in circular 
wells with a radius of r/ = 100 /im and a depth of hf = 
6 /im were investigated. Approximately the same ratio depth 
to radius is typical for the drops of biological fluids used for 
medical tests [7, 8]. 



Moreover, we will suppose that evaporation is a steady state 
process. This assumption is valid e.g. for evaporation of the 
drops of aqueous solutions under room temperature and nor- 
mal atmosphere pressure, i.e. for the typical experimental con- 
ditions. Apex dynamics is rather slow [1] 



H{t) = (7.38/xm- 6.13/xm) - (0.042^m/s)i. 



In the cases of medical tests, the typical velocity of drop apex 
is approximately 1 mm/h. 
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Mass conservation gives the following equation for height 
averaged radial velocity [2] 



(«r(r,t)> = 

1 



prh 



iMh i + 



dh 
Or 



dh 

P ~dl ' Tdr ' ^ 



where h — h(r, t) is the drop shape, p is the density of solution, 
j(r,t) is the evaporation rate defined as the evaporative mass 
loss per unit surface area per unit time. 

Considering the thin droplets only (h(r,t) <C rj) we will 

neglect everywhere h 2 (0, t) and (f^) 2 in compare with 1. We 
will utilize approximative equation for air-liquid interface 



(2) 



h(r,t) = h(0,t) [ 1- ( — 



where r / is radius of circular well. 

Since the drop is thin and the contact angle is small, we will 
use the expression 



Jo 



1 - 



(3) 



for the evaporation rate, which has a reciprocal square-root di- 
vergence near the contact line [6]. This expression can be de- 
rived from Laplace equation (see [2, 6] for details). This func- 
tional form for vapor flux is widely used, in particular, it was 
utilized in paper [1]. 

Some quantities of interest can be expressed analytically 
exploiting (01. Velocity of drop apex decrease is 



dh(0,t) = _4jo 
dt p 
Height averaged radial velocity is 



(v r (x,t)) 



Jo 



(1 - VI - x 2 - 2x 2 



(4) 



(5) 



px(h o + L(0,t)(l-x 2 )) ' 
where x = — , L = — , hn — — . ht is the height of vertical 

re ' r f r f J 

wall of the well. 

Exploiting of Eq.[3]leads to a singularity for both vapor flux 
and gradient of the height averaged velocity of liquid inside a 
drop at the edge of droplet. A smoothing function may be used 
to eliminate physically senseless divergency [9] 

jo 1 — cxp(— mvi — x 2 ) 
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1 — exp(— m) 



(6) 



where m is an adjustable constant. In this case, height averaged 
radial velocity can be written analytically 



(v r (x,t)) = 



Jo 



x(h + L(x,t))m(l~e~ m ) 
2 (l-m-e- m ) i I 



m 



(7) 



where L(x, t) = (l(0, 0) + 4^i=^> j (l - x 2 ) . 
Another evaporative flux functions was proposed in [10] 
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K + L(x, t) 



(8) 



where the constant K measures the degree of non-equilibrium 
at the evaporating interface and is related to material proper- 
ties. K — > corresponds to a highly volatile droplet. The limit 
K — > 00 corresponds to a nonvolatile droplet. Analytical ex- 
pression for height averaged velocity inside circular well can 
be obtained if vapor flux is described by Eq.[8] 



(v r (x,t)) 
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2x(h a + L(0,t)(l- x 2 )) 
( 1 A L(0,t)x 2 



\L(0,t) 



In 1 - 



K 



L(0,t) 



(L(0,t)-2(K + h a )y 
(K + h Q f 



(9) 



where L(0, t) = 2(K+h Q )+(L a -2(K+h )) cxp ( W J ^ ? 

Deegan et al. [2] demonstrated experimentally, that if evap- 
oration is greatest at the center of the droplet, a uniform dis- 
tribution of colloidal particles remained on the substrate. Sim- 
ulations by [3] confirmed, that as fluid is lost from the center 
of the droplet, an inward flow develops to replenish the evapo- 
rated fluid. The evaporative flux function was proposed by Fis- 
cher [3] to mimic evaporation that is concentrated at the center 
of the droplet 
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L(0,t) 



(10) 



where A is an adjustable constant. 

We derived the analytical expression for height averaged 
velocity for evaporation function given by Eq. [10] 



(v r (x,t)) = 

jo (l - e~ Ax " + (1 - e- A ) x 2 (x 2 - 2)) 



2Ax(h + L(x,t))L(0,t) 



, (ID 



where 



L(x,t) 



(l-x 2 )[ X l(ho + Lof-4 Joil - A e ~ A)t -ho 



Fig.|2]demonstrates our calculations for all described above 
evaporative functions. In all figures, the plots corresponds to: 
initial stage (air-liquid interface is convex); air-liquid interface 
is flat; air-liquid interface is concave with the same curvature 
as at the initial stage; 90 % of time, when air-liquid interface 
touch the bottom of the well. 
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Fig. 2. Height averaged radial velocity for different modes of evaporation: a) l[3}; b) {6]l if m — 5; c) l[8}; d) d 10b - For additional explanations 
see the main text. 



2.2 Velocity field 

Velocity fields calculated for sessile droplets from Laplace equa- 
tion [11] and from Navier-Stokes equations [3] are very simi- 
lar. We hope that ideal liquid is reasonable assumption for our 
object of interest. 

Cylindrical coordinates (r, ip, z) will be used, as they are 
most natural for the geometry of interest. The origin is chosen 
in the center of the circular well on the substrate. The coordi- 
nate z is normal to the substrate, and the bottom of the circular 
well is described by z = 0, with z being positive on the droplet 
side of the space. The coordinates (r, tp) are the polar radius 
and the angle, respectively. Due to the axial symmetry of the 
problem and our choice of the coordinates, no quantity depends 
on the angle ip, in particular u = u(r, z). 

The mass flux at the vertical wall of circular well is absent, 
hence 

£ =0, (12) 

or 



where r/ is radius of the circular well. 

Moreover, there is physically obvious relation 



du 
dr 



= 0. 



(13) 



The bottom of the circular well is impermeable, hence 

du 



0. 



(14) 



z=0 



The mass flux inside the droplet near the air-liquid interface 
f(r) is connected with vapor flux 



du 

dz 



= /(r). (15) 

free surface 

Laplace equation in cylindrical coordinates for the axially 
symmetric case is written as 

Id/ du\ d 2 "' 



r dr \ dr 



dz'< 



= 0. 



(16) 
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Boundary problem (|12t . (|14t . for equation (TToT i 

u(r, z) = } a n cosh z Jo r 

±1 V r f J \ r f J 



(17) 



where J m (r) is the Bessel function of the first kind, order m, 

fin are the real zeros of Bessel Ji(r) function: J\{r) = 0. 
Note that condition ( fl~3T > is satisfied automatically. 

Taking into account relation between velocity and potential 
v = — grad it, we can find radial component of velocity 

(J£) \ M (i) (M) \ 
(v r (r,z)) = y^a„cosh z -^-Ji -^—r . 

\ r f J r f \ r f J 

Height averaged radial velocity, i.e. the same velocity ex- 
amined in [1], can be written as 



Mr)) 



v r (r, z) dz 



— y~] a n sinh -^—ht J x -^—r 



Coefficients a n can be obtained from (fT~5T >. For simplicity 
we will assume that air-liquid interface is flat. 



du 

dz 



a n — amh — h f J — r 



f(r). (18) 



We computed velocity fields for all evaporation laws de- 
scribed in Sec. 12. 11 The results show that current inside circular 
well is horizontal excluding rather narrow region near the wall 
and in the center of the well for the case of practical interest 
(/i/r~0.1). 

Many authors supposed that evaporation induced flow in- 
side the droplets is independent of its height. Our simulations 
confirmed validity of this very wide used approach for thin 
droplets. 

The authors are grateful to the Russian Foundation for Basic Research 
for funding this work under Grant No. 06-02- 16027-a. 
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We introduce notation 



f(r)r dr — F(r). 



V* c n Ji -^—r 
±1 V r f j 



F(r) 



is a Fourier-Bessel series of F(r) / r, where 
c n = a n sinh h f . 



(19) 



(20) 



Hence, 



/u)\ 7 ET\ I Jl (— r ) F ^ dr - (21) 

We need to determine f(r) to find the unknown coeffi- 
cients. 

Mass conservation leads to the relation 



du 

dz 



dL{r,t) _ J(r,t) 



(22) 
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Fig. 3. Velocity field inside circular well for differen evaporative modes. The air-liquid interface is supposed to be flat. From top to bottom {3), 
(|6j with m = 5, d 1 Ob - Left column corresponds to ratio hf jrj = 0.0618, right column corresponds to ratio hf jrf — 0.99 



